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Abstract 

We present a calculation of the charm and gluon fragmentation contributions to 
inclusive J/ip and ip' production at large transverse momentum at the Tevatron. For ip 
production, we include both fragmentation directly into ip and fragmentation into Xc 
followed by the radiative decay Xc ~~ ^ ip+'J- We find that fragmentation overwhelms the 
leading-order mechanisms for prompt ip production at large pr, and that the dominant 
contributions come from fragmentation into Xc- Our results are consistent with recent 
data on ip production from the CDF and DO experiments. In the case of prompt ip' 
production, the dominant mechanism at large pt is charm fragmentation into ip' . We 
find serious disagreement between our theoretical predictions and recent ip' data from 
the Tevatron. 

1 on leave from the Dept. of Physics and Astronomy, Northwestern University, Evanston, IL 60208 



1 Introduction 

The study of charmonium production in high energy hadronic collisions provides an impor- 
tant testing ground for perturbative quantum chromodynamics (QCD). The J/tp and ip' 
states are of particular interest since they are produced in abundance and are relatively easy 
to detect at a collider such as the Tevatron. In earlier calculations of direct charmonium 
production at large transverse momentum (jpr) in pp collisions |1[ , it was assumed that the 
leading-order diagrams give the dominant contributions to the cross section. These calcu- 
lations did not reproduce all aspects of the available data @, §, f§], suggesting that there 
are other important production mechanisms. It was pointed out by Braaten and Yuan || 
in 1993 that fragmentation processes, while formally of higher order in the strong coupling 
constant a s , will dominate at sufficiently large pt- Explicit calculations of the contribution 
to ip production at the Tevatron from the fragmentation of gluons and charm quarks revealed 
that fragmentation dominates over the leading-order gluon-gluon fusion mechanism for p? 
greater than about 6 GeV ||. In addition to being directly produced, the ip signal is also 
fed by the \c states through the radiative decay \c — > i> + 7- In this paper, we present more 
complete calculations for ip production, including the effects of Xc's that are produced by 
fragmentation, and compare the results with recent data on prompt and inclusive ip produc- 
tion from the CDF and DO experiments. We also compare our calculations of ip' production 
with recent data from CDF. 

While this paper was being written, similar work on -^-production was presented in 
a paper by Cacciari and Greco [0. 

2 Fragmentation Formalism 

Factorization theorems of perturbative QCD indicate that the inclusive production of a 
hadron at large Pt is dominated by fragmentation. Fragmentation is the production of a 
parton with large transverse momentum which subsequently forms a jet containing the de- 
sired hadron. In the case of pp collisions, the fragmentation contribution to the cross section 
can be expressed as a convolution of parton distribution functions, hard-scattering cross 
sections, and fragmentation functions. Taking ip production to be specific, the differential 
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cross section can be written as: 

da(pp -> ip(p T , y) + X) = J2 [ dz d(J {PP -> «(— , V) + x , /^frag 2 ) D^(z, /i frag 2 ), (1) 

where z is the longitudinal momentum fraction of the if) relative to parton i, y is the ra- 
pidity of the if), and Di->^(z, /if rag 2 ) is the fragmentation function. The dependence on the 
fragmentation scale /if rag cancels between the two factors only after inclusion of all orders 
of the perturbative expansion. The differential cross section on the right side of ([!]) can in 
turn be written as the convolution of parton distributions fj/ p and fk/p in the proton and 
antiproton with hard-scattering differential cross sections da for the parton subprocesses 
j + k — > i + X. In low-order calculations, the fragmentation scale /Xf rag , the factorization 
scale /if , and the renormalization scale /ir, should all be chosen on the order of Pt/z, the 
transverse momentum of the fragmenting parton. The dominant contributions to Eq. ([!]) 
come from gluon fragmentation and charm fragmentation. Fragmentation of light quarks 
will only contribute at low z, and is therefore strongly suppressed by the rapidly falling 
spectrum of the final-state partons from the hard-scattering process. The inclusive ip signal 
is also fed by production of the P-wave Xc states followed by the radiative decay \c — > i> + 7- 
The fragmentation contributions can be obtained by multiplying the Xc production cross 
sections analogous to ([!]) by the radiative branching fractions 0.007, 0.27, and 0.135 for x c0 , 
Xd, and Xc2, respectively. 

The fragmentation functions D(z) for charmonium production can be calculated 
within perturbative QCD 0, The relevant fragmentation functions for the production of 
the S-wave and P-wave states have all been calculated to leading order in a s . In this paper 
we use the fragmentation functions for g — * if) ||, c — > if) [0, g — ► Xc [pT| , c — > Xc |]T2 



and 7 — > if) |^3| . The perturbative calculations give the fragmentation functions at an initial 
scale (jLq of order m c . We take this initial scale to be /io = 2m c for gluon and photon frag- 
mentation, and = 3m c for charm fragmentation. The fragmentation functions are then 
evolved up to the scale /if rag = 0(px/z) set by the transverse momentum of the fragmenting 
parton using the Altarelli-Parisi evolution equations: 

^D^(z,^) = g£ J^PMv) iWv, A (2) 



We included in Eq. (|2]) only the P gg splitting term for gluon fragmentation and only the P cc 
term for charm fragmentation. With the exception of D g ^, the inclusion of the off-diagonal 
Altarelli-Parisi kernels would give corrections that are safely negligible within the accuracy 
of a leading-order (LO) calculation. In the case of D g ^, the P gc term in (H) is important 
because the initial fragmentation function D c ^(z, /ig) is more than an order of magnitude 
larger than D g ^(z, fj^) ||. In the present analysis, we take this effect into account by 
including the next-to-leading-order (NLO) corrections to the hard-scattering cross sections 
da for inclusive charm production |TJ]]. These higher order terms include the numerically 
important contributions from final state gluon splitting into cc pairs. Having included this 
effect in the hard-scattering cross sections, we do not have to account for them via the g — > cc 
Altarelli-Parisi kernel as in Ref. 0. Higher-order terms from off-diagonal Altarelli-Parisi 
evolution, such as g — > cc splitting from secondary gluons in the gluon-jet cascade, contribute 
to the overall charmonium multiplicity inside the jet, but the enhancement comes only in 
the small z region and it can therefore be safely neglected in the LO px distribution. 

We now discuss the uncertainties in the calculations of the fragmentation contribu- 
tion to the differential cross section for charmonium production. We first consider the initial 
fragmentation functions D(z,/1q). They have been calculated only to leading order in a s 
and to leading order in a nonrelativistic expansion. Based on the NLO calculations of the 
annihilation rates for charmonium, we anticipate that NLO corrections to the fragmentation 
functions may be as large as 50%. A rough estimate of the size of relativistic corrections is 
the average value of v 2 /c 2 in potential models, which is about 30%. The normalizations of 
the initial fragmentation functions are calculated in terms of a s , m c , and various nonpertur- 
bative matrix elements. For the charm quark mass, we use the value m c = 1.5 GeV. While 
m c appears in the fragmentation functions raised to the third or fifth power, this is not a 
large source of uncertainty since roughly the same power appears in the quantities that are 
used to determine the nonperturbative matrix elements. The S-wave fragmentation functions 
depend on the wavefunction at the origin, which we take to be 1^(0) | 2 = °- 7 GeV3 - This 
value is obtained from the electronic width of the ip, including the effect of the NLO pertur- 
bative correction, which is about 50%, but not taking into account relativistic corrections. 
The P-wave fragmentation functions depend on two nonperturbative parameters. For the 
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derivative of the radial wavefunction at the origin, we use the value |i?^. c (0)| 2 = 0.053 GeV 5 . 
This value is determined from the annihilation rates of the \ c states , neglecting the as- 
yet-uncalculated NLO perturbative corrections as well as relativistic corrections. The least 
well-determined parameter in the P-wave fragmentation functions is a parameter H' s associ- 
ated with the color-octet mechanism for P-wave production This parameter is poorly 
constrained, lying in the range 2.2 < H' 8 < 25 MeV [17|. We use the value H' B = 3 MeV, 



which is consistent with measured branching fractions for B mesons into the \ c states [16 



We next consider errors due to the evolution of the fragmentation functions from the 
initial scale fio up to the scale Pt/z. The Altarelli-Parisi equations Eq. (^) break down in 
the small-;? region, due to large logarithms of 1/z in the perturbation expansion. The most 
dramatic effect of this breakdown is an unphysical divergence in the gluon multiplicity at 



low momentum fraction . This leads to a correspondingly large overestimate of the gluon 
fragmentation functions into charmonium at small values of z. We do not expect this to be 
a serious problem at the values of px considered in this paper, because the cross sections are 
dominated by larger values of z due to the steeply falling spectrum of the hard partons. We 
find empirically that the Monte Carlo calculations presented below rarely sampled values of 
z smaller than 0.1. 

While the divergence of the gluon multiplicity at small z may not in itself be a prob- 
lem, it is a symptom of a deficiency of Altarelli-Parisi evolution that also has consequences at 
larger values of z. In particular, the naive Altarelli-Parisi equations do not respect the phase- 
space constraint D g ^(z, fi 2 ) = for z < M 2 /fi 2 . The implementation of this constraint 
would slow down the evolution of the fragmentation function by delaying the depletion of the 
large- z fragmentation region. Since the spectrum of gluons and charm quarks falls rapidly 
with Pt/z, a proper treatment of the large- z region can have a significant effect on the cross 
section. A more accurate treatment would be to use the following system of equations for 
the evolution of the gluon fragmentation functions: 

* 2 da* rl 



D(z^ 2 ) = f %/ 1 ^C(y,g 2 ;^)d(z/y,g 2 ) (3) 
Jm% q 2 Jz y 



^tt^^'^^o) = it- I — p ag(y) G ( z /y^ 2 ;yvl) (4) 

oyi z Z7T Jz y 



where G(z, /x 2 ; /i 2 ,) is the distribution of gluons of virtuality /x inside a gluon of virtuality /xo, 
which is subject to the boundary condition G(z, /X 2 ,; /X 2 ,) = 8(1— z). The function d(z, q 2 ) || is 
the LO probability that a gluon of virtuality g decays to a ^ carrying longitudinal momentum 
fraction z in the infinite momentum frame. The initial fragmentation functions in our present 
treatment are given by D(z, /Xq) = fo°(ds/s)d(z, s). The system of equations (|3|) and (^) 
gives the correct ip multiplicity inside a gluon jet | I9| , after inclusion of small- 2 coherence 



effects |I8|. It is straightforward to show that this system is equivalent to the following 



nonhomogeneous evolution equation: 

^ 2 ^- 2 D(z^ 2 ) = d(z^ 2 ) + ^( ld JL Pgg (y) D(z/y, W 2 ), (5) 

together with the boundary condition D(y,fi 2 = M 2 ) = 0. This evolution equation respects 
the phase space constraint, as can be easily checked |IS|. A thorough study of this general- 



ized evolution equation and its consequences for charmonium production will be presented 
elsewhere. 

The evolution equation presented above also solves a problem involving threshold 
effects in our fragmentation functions. The px values for which the fragmentation contri- 
butions become important are, in fact, too small to be considered in the asymptotic regime 
where threshold effects can be ignored. In the calculations of the initial fragmentation func- 
tions _D(z,/Xq), the assumption /x 2 , S> 4m 2 /z was used to obtain simple analytic results. This 
assumption is not really compatible with the subsequent identification /x = 2m c . The re- 
sulting error may be negligible after evolution to asymptotically large /x, but it is probably 
significant at the scales needed in our calculations. The evolution equation (f|) treats thresh- 
old effects consistently. In the present analysis, we estimate the error due to our treatment of 
threshold effects by determining the change in the cross section that results from increasing 
the value of /xo by a factor of 2. 

There are several other sources of uncertainty. The uncertainties due to the choice 
of the renormalization scale, the factorization scale, and the fragmentation scale can be 
estimated by varying these scales by factors of 2. The error from the parton distributions 
can be estimated by repeating the calculations using different parton distributions. Another 
source of error at small pt is the neglect of the intrinsic transverse momentum of the partons 



in the proton and antiproton [EH, |21[ . The effect of the intrinsic transverse momentum is most 



significant for partons with very small longitudinal momentum fraction x. It should not be 
important at large pt because ip production at large pr does not probe deeply into the small- 
x region of the parton distribution. Finally there is the NLO perturbative correction to the 
hard-scattering cross section da, which we only included in the case of charm fragmentation. 
The NLO corrections to the gluon p? spectrum have been evaluated in ref. , where it 
was shown that they increase the rate for x production from fragmentation by about 50%. 
Needless to say, a full NLO calculation of both the hard-scattering cross sections and the 
fragmentation functions would be highly nontrivial. 

Considering all the uncertainties discussed above, we believe that the error in our 
fragmentation calculations can easily be larger than a factor of 2, but it is definitely less than 
an order of magnitude. One should of course keep in mind that in addition to fragmentation, 
which must dominate at sufficiently large p?, there are other contributions suppressed by 
factors of m 2 c /p T that may be important at the values of pt that are available experimentally. 



3 Results and discussion 

In figure 1 we plot the individual contributions to the differential cross section for prompt ip 

production as a function of pr- We include results for both the fragmentation contributions 

and the leading-order contributions^. Note that each of the \ c production curves is a sum 

over the three P-wave states Xco, Xa, and \c2- We used the MRSDO parton distribution set, 

and chose the renormalization (/xr), factorization (/ip) and fragmentation (/if rag ) scales to 

be same, and equal to the transverse momentum of the fragmenting parton, Pt=Pt/ z - In 

order to compare with available data from the Tevatron, we imposed a pseudorapidity cut 

of 1 7/| < 0.6 on the ip. It is evident from the graph that fragmentation dominates over the 

leading-order mechanisms for all values of pt for which the fragmentation approximation is 

reasonable, namely for p? greater than about 5 GeV. The dominant production mechanism by 

an order of magnitude is gluon fragmentation into \ c followed by its decay into ip. Note that, 

aside from photon fragmentation which is dominated by quark-gluon initial states, all the 

2 The results for direct i/j inclusive pt distributions given in ref. [Q are incorrect, due to a coding error in 
the choice of [1r and ^f- 
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fragmentation contributions have the same pr-dependence. The leading-order Xc production 
falls off more rapidly with px, and that leading-order ip production falls off still more rapidly. 
This pattern simply reflects the pr-dependence of the underlying hard-scattering processes: 
da/dp T scales like l/p T for fragmentation, \jp% for leading-order Xc production, and l/p T 
for leading-order ip production. 

In figure 2, the sum of the fragmentation contributions (two solid curves) and the sum 
of the leading-order contributions (two dashed curves) are compared with preliminary CDF 
data for prompt ip production ||22|| . The contribution to ip production from 6-hadron decays 
has been removed from the data via detection of the secondary vertex from which the ^'s 



originate |[22|| . The upper and lower curves in figure 3 were obtained by varying the scales 
/xr, /x f and /if rag used in the calculation, in order to provide an estimate of the systematic un- 
certainty associated with the LO calculation. The upper curve corresponds to /xr=/xf=-Pt/2 
and /Xf rag =max(P T /2,/i ), while the lower curve is obtained for /XR=/x F =/Xf rag =2P T . The 
cross-over of the curves at small pt is due to the rapid growth of the parton distribution 
functions with increasing scale, and should be considered an artificial reduction of scale sen- 
sitivity. The errors from varying the parton distributions are only about 10%. We estimate 
the error from threshold effects in the fragmentation functions by increasing /xo by a factor 
of 2. The effect of this is a decrease of the fragmentation contribution, by a factor of 2, at 
Pt = 3 GeV, and an increase, by a factor of 2, at px = 20 GeV. Another large uncertainty 
comes from the color-octet parameter H f 8 in the fragmentation functions for g — > Xc- Chang- 
ing H f 8 by a factor of two changes our results by a factor of 2 at the largest px available. 
While the shapes of the leading-order curve and the fragmentation curve are both consistent 
with the data over the range of px that is available, the normalization of the leading-order 
contribution is too small by more than an order of magnitude. The fragmentation contribu- 
tion has the correct normalization to within a factor of 2 or 3, which can be easily accounted 
for by the uncertainties discussed above. We conclude that the fragmentation calculation is 
not inconsistent with the CDF data on prompt ip production. 

We also present in fig. 3 a comparison of the theoretical predictions for fully inclusive 
i/j production (including those from 6-hadron decays) with data from CDF |^2[ and DO |23 |. 
The theoretical contribution from 6-decays used here was evaluated at NLO as in Ref. H. 
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The predictions are consistent with the data. 

We next consider the production of if)', which should not receive any contributions 
from higher charmonium states. The if)' fragmentation contribution can be obtained from 
the g — *> if), c — > if), and 7 — > if> fragmentation contribution simply by multiplying by the 
ratio of the electronic widths of the if>' and if>. The total fragmentation contribution (two 
solid curves) and the leading-order contribution (two dashed curves) are shown in figure 4, 
along with the preliminary CDF data p4 |. Again the contribution from 6-hadron decays has 
been subtracted using the secondary vertex information. The pairs of curves correspond to 
the same choices of scales as in figure 2. The dominant production mechanisms are gluon- 
gluon fusion for px below about 5 GeV, and charm quark fragmentation into if)' for larger 
Pt- The leading-order curve falls much too rapidly with p? to explain the data, but the 
fragmentation curve has the correct shape. However, in striking contrast to the case of if) 
production, the normalization of the fragmentation contribution to if)' production is too small 
by more than an order of magnitude. That there is such a large discrepancy between theory 
and experiment in the case of if)', but not for if), is extremely interesting. It suggests that 
there are other important mechanisms for production of S-wave states at large px beyond 
those that have presently been calculated. While such processes would certainly affect if) 
production as well, their effect may not be as dramatic because of the large contribution 
from x c -production in the case of the if). 

One possible such mechanism is the process gg — > ipgg, with a gluon exchanged in 
the t-channel. This is a subset of the NLO corrections to the process gg — > if>g for which the 
hard-scattering cross section has a pr-dependence that is intermediate between the leading- 
order diagrams and the fragmentation contribution. It is easy to verify that da/dp T for this 
process scales asymptotically like ot s A /pr 6 , compared to a s 3 /pT 8 for gg — > ipg and compared 
to a s 5 /pT A for the NNLO correction, which includes gluon fragmentation. The shape of the 
Pt distribution from this process should be similar to that for leading-order \ c production, 
which is compatible with the data. Whether the normalization agrees can only be determined 
by explicit calculation. Such a calculation is in progress. 

Note that, like the decay of a 6-hadron, the fragmentation mechanism produces if)'s 
inside a jet of light hadrons. A non-isolation cut on the if) can therefore not be used to 
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tag ^'s coming from B meson decay. One might hope that an isolation cut could be used 
to separate prompt ^'s produced by the leading-order mechanisms from those produced by 
fragmentation. However, we found in our calculations that the average value of z for a ip 
produced by fragmentation is (z) = 0.7. This means that the remaining partons in the jet 
containing the ip share on average less than 1/2 the energy of the tp, and thus may often be 
too soft to be detected reliably. 

4 Conclusion 

We have calculated the cross section for production of prompt ip at large p T at the Tevatron 
including the effects of gluon and charm fragmentation. The largest contributions by an order 
of magnitude come from gluon fragmentation into Xc followed by the decay Xc ip + 7- 
The results of this calculation are not inconsistent with preliminary CDF data, given the 
uncertainties in the fragmentation functions. In contrast, the leading-order mechanisms give 
a result that is nearly an order of magnitude too small. 

We have also calculated inclusive ip' production at large px at the Tevatron. We find 
that fragmentation dominates for p T greater than about 6 GeV, with the largest contribu- 
tion coming from charm fragmentation into ip'. Comparing the results of the calculation to 
preliminary CDF data, we find that the cross section is too small by more than an order of 
magnitude, even after including the fragmentation mechanism. Thus while the fragmenta- 
tion mechanism may provide an explanation for the observed rate of prompt ip production 
at large pt, it does not seem to explain the existing data for ip' . We discussed the possibility 
of additional mechanisms for charmonium production at large px beyond those that have 
presently been calculated. Furthermore, we discussed some of the intrinsic theoretical uncer- 
tainties present in these calculations and anticipated some studies that can be undertaken 
to improve the present theoretical framework. 
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Figure Captions 



1. Contributions to the differential cross section for inclusive ip production at the Teva- 
tron: fragmentation into ip (solid curves), and the leading order contributions (dashed 
curves) . 

2. Preliminary CDF data for prompt ip production (O) compared with theoretical predic- 
tions of the total fragmentation contribution (solid curves) and the total leading-order 
contribution (dashed curves). 

3. Total ip production: CDF (O) and DO (O) data compared to theoretical curves for 
prompt ip production (solid curves), and theoretical predictions for 6-hadron decays 
(dashed curves). 

4. Preliminary CDF data for prompt ip' production (O) compared with theoretical predic- 
tions of the total fragmentation contribution (solid curves) and the total leading-order 
contribution (dashed curves). 
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